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ABSTRACT 

We explore the gravitational instability of clumpy and turbulent gas discs, taking into 
account the Larson-type scaling laws observed in giant molecular clouds (GMCs) and 
Hi, as well as more general scaling relations. This degree of freedom is of special 
interest in view of the coming high-z ISM surveys, and is thus potentially important 
for understanding the dynamical effects of turbulence at all epochs of galaxy evolution. 
Our analysis shows that turbulence has a deep impact on the gravitational instability of 
the disc. It excites a rich variety of stability regimes, several of which have no classical 
counterpart. Among other diagnostics, we provide two useful tools for observers and 
simulators: (1) the stability map of turbulence, which illustrates our stability scenario 
and relates it to the phenomenology of interstellar turbulence: GMC/Hi observations, 
simulations and models; (2) a Toomre-like stability criterion, Q > Q, which applies to 
a large class of clumpy/turbulent discs. We make specific predictions about GMC and 
cold-Hl turbulence, and point out the implications of our analysis for high-z galaxy 
surveys. 
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1 INTRODUCTION 

Toomre's (1964) stability criterion, Q > 1, is one of the pil- 
lars of disc galaxy dynamics (see, e.g., Binney & Tremaine 
2008). It is used in a wide variety of contexts; in star forma- 
tion for example, where the gravitational instability of the 
interstellar gas plays a critical role (e.g., Quirk 1972; Ken- 
nicutt 1989; Martin & Kennicutt 2001; Schaye 2004, 2008; 
Burkert 2009; Elmegreen 2009). 

So why introduce a new stability criterion? Because be- 
hind Toomre's criterion there is one fundamental assump- 
tion: the medium is approximately in equilibrium, with well- 
defined surface density S and velocity dispersion u. But this 
is far from being true in the clumpy and turbulent interstel- 
lar gas, where such quantities depend strongly on £, the size 
of the region over which they are measured. In fact, a fun- 
damental aspect of interstellar turbulence is the existence of 
density- and velocity-size scaling laws: E oc and a oc £^ 
(see, e.g., Elmegreen & Scalo 2004; McKee & Ostriker 2007). 
In giant molecular clouds (GMCs), the scaling exponents are 
a « and b ~ 5 (e.g., Larson 1981; Solomon et al. 1987; 
Bolatto et al. 2008; Heyer et al. 2009; Hughes et al. 2010). 
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In the H I component, density and velocity fluctuations seem 
to have a Kolmogorov spectrum up to galactic scales: a ~ | 
for ^<10 kpc, and & ~ i for ^<1 kpc (e.g., Lazarian & 
Pogosyan 2000; Elmegreen et al. 2001; Begum et al. 2006; 
Kim et al. 2007; Dutta et al. 2008; Roy et al. 2008; Dutta et 
al. 2009a, b). Note, however, that the uncertainties are large, 
especially in the Hi case. Further evidence for ISM turbu- 
lence is provided by the 'Big Power Law in the Sky', i.e. the 
fact that electron density fluctuations show a Kolmogorov 
spectrum over a wide range of scales: from 10"^ km up to 10 
pc (Armstrong et al. 1995; see Chepurnov & Lazarian 2010 
for the most recent determination of the upper scale). 

Clumpy and turbulent gas is also observed in high- 
redshift star-forming galaxies, where it dominates the mor- 
phology and dynamics of the disc (e.g., Wadadekar et al. 
2006; Elmegreen et al. 2007; Genzel et al. 2008; Shapiro 
et al. 2008; Forster Schreiber et al. 2009). Coming surveys 
will tell us how E and a scale with £ at high z, and thus 
how turbulence develops in disc galaxies. This is one of the 
hot topics in modern astrophysics (e.g., Krumholz & Burk- 
ert 2010), and several state-of-the-art simulations have al- 
ready been designed for such a purpose (e.g., Wada et al. 
2002; Kim & Ostriker 2007; Levine et al. 2008; Agertz et 
al. 2009a, b; Bournaud & Elmegreen 2009; Dekel et al. 2009; 
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Tasker & Tan 2009). In order to interpret such data cor- 
rectly, one must understand in detail how turbulence affects 
the (in)stability of the disc. 

These facts all together motivate a thorough investiga- 
tion of the problem. 

Elmegreen (1996) assumed Larson- type scaling rela- 
tions, E oc and a oc and investigated the case o = — 1 
and 6 = i. He found that the disc is always stable at large 
scales and unstable at small scales. To the best of our knowl- 
edge, this was the only theoretical work devoted to the grav- 
itational instability of clumpy/turbulent discs before ours. 
In contrast, several investigations focused on the effects of 
turbulence on Jeans instability (e.g., Bonazzola ct al. 1987; 
Just et al. 1994; Vazqucz-Scmadcni & Gazol 1995). 

The goal of our paper is to explore the gravitational 
instability of clumpy /turbulent discs, spanning the whole 
range of values for a and b. Among other diagnostics, we 
provide two useful tools for observers and simulators: 

(i) The stability map of turbulence. Using a and b as coor- 
dinates, we illustrate the natural variety of stability regimes 
possessed by such discs, and populate this diagram with ob- 
servations, simulations and models of interstellar turbulence 
(Fig. 1). 

(ii) A Toomre-like stability criterion. We show that in our 
map there is a densely populated domain where the stability 
criterion is of the form Q > Q, and determine the stability 
threshold Q as a function of a, b and of the scale at which Q 
is measured [Eq. (17) and Fig. 2]. If our criterion is fulfilled, 
then the disc is stable at all scales (the case investigated by 
Elmegreen 1996 lies in another stability regime). 

The rest of the paper is organized as follows. In Sect. 
2, we determine a dispersion relation that takes fully into 
account the scaling laws of interstellar turbulence as well as 
the thickness of the gas layer. We then surf through the var- 
ious stability regimes. Our Toomre-like stability criterion is 
given in Sect. 2.7, together with other important diagnos- 
tics: the most unstable scale and its growth rate. In Sect. 3, 
we discuss the stability map of turbulence. We then make 
specific predictions about strongly and weakly supersonic 
turbulence. In Sect. 4, we explore the densely populated 
Toomre-like domain and illustrate the stability diagnostics 
in a number of cases, taking into account the saturation of 
density and velocity at large scales. In Sect. 5, we draw the 
conclusions. 



2 STABLE OR UNSTABLE ? 

2.1 The turbulent dispersion relation 

Consider a gas disc of scale height h, and perturb it with 

axisymmetric waves of frequency and wavenumber k. The 
response of the disc is described by the dispersion relation 

2 2 SttGE k 



1 + kh 



1 2,2 



(1) 



where k is the epicyclic frequency, E is the surface density 
at equilibrium, and a is the sound speed (Romeo 1990, 1992, 
1994; see also Vandervoort 1970). So the three terms on the 
right side of Eq. (1) represent the contributions of rotation, 
self-gravity and pressure. For fcft. -C 1, Eq. (1) reduces to 
the usual dispersion relation for an infinitesimally thin gas 



disc (see, e.g., Binney & Tremaine 2008). For fc/i S> 1, one 
recovers the case of Jeans instability with rotation, since 
E//i « 2p. In other words, scales comparable to h mark 
the transition from 2D to 3D stability. Note that we can 
encapsulate the effect of thickness in a single quantity and 
rewrite Eq. (1) as 



— 27rGEefi k + a'^ k^ . 



(2) 



where Eeff = E/(l -|- kh) is the effective surface density. 
Eeff and a are two important quantities, which we discuss 
in detail below. 

The effective surface density. What is the relation between 
Ecfi and p? Since Eeff = E/(l + kh) and E ~ ph, we find 
that Eefi; ~ ph if fc/i<l, and Eefi ~ p/k otherwise. This 
means that the observational counterpart of Eeg is the mass 
column density measured over a region of size £ = 1/k. In 
the cold ISM, which is highly supersonic and hence strongly 
compressible, the amplitude of density fluctuations is typ- 
ically much larger than the mean density. Density fluctu- 
ations have a power-law energy spectrum, Ep{k) oc fe"'", 
which means p oc k~^'^~^^^^ (see, e.g., Elmegreen & Scalo 
2004). A power-law spectrum is then imprinted on Eefi: 



Eeff 



(3) 



Using the relation between Eefi and p that we have found 
above, we can relate a to r: o = | (r — 1) for kh < 1, while 
a = \{r + 1) for fe/i> 1. In the warm ISM, which is tran- 
sonic or subsonic and hence weakly compressible, the den- 
sity contrast is typically small so Eeff is dominated by the 
mean density, as in the limiting case of a non-turbulent disc: 
Eefi = constant for small k, while Eeff oc fc^^ for large k. 
Hereafter we will omit the subscript 'eff', unless otherwise 
stated. 

The sound speed. The observational counterpart of a is 
the ID velocity dispersion measured over a region of size 
i = 1/k: = (Tjijer + crtur(fe), where uther and utur are the 
thermal and turbulent ID velocity dispersions. Velocity fluc- 
tuations have a power-law energy spectrum, Ev{k) oc fc"", 
which means atm oc fc"'-""^'''^ (see, e.g., Elmegreen & Scalo 
2004). Both thermal and turbulent motions tend to support 
the gas against gravitational instability. However, as pointed 
out by the referee, turbulent support results in relations that 
are true only in a statistical sense, with individual clouds col- 
lapsing even when the statistical criterion for stability is well 
satisfied. In the cold ISM, a has a power-law dependence on 
k since it is dominated by cTtur(/c): 



(4) 



where b — i (s — 1) and s is the velocity spectral index. In 
the warm ISM, a{k) deviates significantly from a power law 
since fXther is no longer negligible. 

Hereafter wc will only consider the cold ISM (H2 and 
cold Hi), since in such a case we can explicitely take into 
account the power-law scaling of interstellar turbulence via 
Eqs (3) and (4). We will regard the quantity £0 = l/fco 
introduced in those equations as the fiducial scale at which 
the mass column density and the ID velocity dispersion are 
observed. This is also the scale at which Q and other stability 
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quantities are measured. Substituting Eqs (3) and (4) into 
Eq. (2), we obtain the final dispersion relation: 



2u2b ,2(1-6) 



(5) 



Our phenomenological approach diflters significantly 
from the traditional way to include turbulent effects in the 
dispersion relation, which is to identify a with the typical 
ID velocity dispersion observed at galactic scales. Bonaz- 
zola et al. (1987), Just et al. (1994), Vazquez-Semadeni & 
Gazol (1995) and Elmegreen (1996) adopted an approach 
similar to ours for investigating the gravitational instability 
of clumpy /turbulent media. 

Eq. (5) is the starting-point of our stability analysis. As 
in the usual case, stability at all scales requires that o;^ > 
for all fcQ Whether this requirement can be satisfied or not 
depends on the self-gravity and pressure terms, which now 
scale as k^ and k^: 

A = 1 - a , (6) 
5 = 2(1-6). (7) 
In the following, we explore all the various cases. 



2.2 Case a = 1, fe / 1 

Here the self-gravity term is independent of fc (A = 0), like 
the rotation term: 



2 2 o /-fv^ 7 I 2 , 2b , 2(1-6) 



(8) 



c 



For h < 1, uj^ increases with k and tends to C as A; 0. 
Hence the sign of C determines whether the disc is stable at 
all scales or not. Stability requires that C > 0. For 6 > 1, 
decreases with k and tends to C as — i- oo. So, again, 
stability requires that C > 0. The stability criterion is then: 



STABLE Vfc 



fco < 



27rGEo 



(9) 



Eq. (9) resembles Toomre's stability criterion for cold discs, 
k < kT, but there is one important difference: Eq. (9) is a 
condition that must be fulfilled by fco to ensure stability for 
all fc, whereas the cold Toomre criterion ensures stability for 
small fc. 



If a < 1, then lo converges to k as fc — >■ and the sign 
of doj^/dfc is equal to the sign of C (for C < 0, cj^ diverges 
as fc — > oo). Hence stability requires that C > 0. If a > 1, 
as fc — >■ cx) and the sign of da;^/dfc 

,2 



then uj'^ converges to 



is opposite to the sign of C (for C < 0, oj diverges as fc — >■ 

0). So, again, stability requires that C > 0. The stability 
criterion is then: 

27rGEo 



STABLE Vfc • 



fco > fcj = 



(11) 



The resemblance between Eq. (11) and the 2D Jeans crite- 
rion, fc > fcj, is superficial. In fact, Eq. (11) ensures that the 
disc is stable at all, rather than small, scales. An analogous 
fact was noted in the context of our first stability criterion 
[Eq. (9)]. 

2.4 Case a = 1, 6=1 

Although this case seems the intersection of Case 2.2 and 
Case 2.3, the stability criterion is not fcj < fco < fci. In fact, 
this case is singular. The self-gravity and pressure terms are 
independent of fc {A — B = 0), like the rotation term, so 
ix)^(fc) is constant: 

uj^ = ~ 27rGEofco -|- ctq^o • 



(12) 



C 



As C is quadratic in fco, the inequality G > can be easily 
solved and the resulting stability criterion is: 



STABLE Vfc 



fco < fc- or fco 
< fco < oo 



> fc4 



if Q < 1 , 
else . 



(13) 



Here fc- and fc+ are related to the Toomre wavenumbcr fci 
(and to the Jeans wavenumber fcj = fcx 4/Q^): 

2 



fc± 



Q 



1-0' 



and Q is the Toomre parameter: 



Q 



Kao 
ttGEo 



(14) 



(15) 



The Q < 1 case of Eq. (13) resembles the corresponding 
Toomre stability condition, but see the remarks following 
Eqs (9) and (11). In contrast, the Q > 1 case is identical to 
Toomre's criterion. 



2.3 Case a / 1, 6 = i (1 + a) 

Here the pressure term has the same fc-dependence as the 
self-gravity term {B = A): 

, l-a 



uj^ = + ((TofcQ^° — 27rGEofco) k^ 



(10) 



G 



2.5 Case a < 1, 6 > i (H- a) 

Now j4 > and B < A. Hence the self-gravity term gets 
dominant for large fc and makes uj^ negative. For small fc, 
uj^ is positive since it is dominated by the pressure term 
{B < 0) and/or the rotation term {B > 0). As the zero(s) 
of LO^{k) can only be determined numerically, case by case, 
we do not give a stability criterion but note that the disc is 
unstable at small scales, like a cold non-turbulent disc. 



^ The dispersion relation assumes that kR 2> 1, where R is the 
radial coordinate (see, e.g., Binney & Tremaine 2008). This lo- 
cal condition is more restrictive than the natural requirement 
fc>27r/L, where L is the size of the gas disc. Since the condi- 
tion above cannot be rigorously included in the stability analysis, 
the usual procedure is to consider all k and to check a posteriori 
whether kR 2> 1 or not. 



2.6 Case a > 1, 6 < i (1 a) 

In contrast to the previous case, A < and B > A. So itP' is 
dominated by the self-gravity term and is negative for small 
fc, while it is dominated by the pressure/rotation term and is 
positive for large fc. Thus the disc is unstable at large scales, 
like a non-rotating non-turbulent sheet. 
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2.7 Case a < 1, h <\{l + a) 



3 THE STABILITY MAP OF TURBULENCE 



When Q < A < B, the response of the disc is driven by pres- 
sure at small scales and by rotation at large scales, while 
self-gravity acts more strongly at intermediate scales. There- 
fore this is a Toomre-like case: w^(fe) has a minimum, which 
determines whether the disc is stable or not. More precisely, 
the minimum of uj^{k) provides three useful pieces of infor- 
mation: the stability threshold, the most unstable scale and 
its growth rate. Such quantities are introduced below, while 
illustrative cases are discussed in Sect. 4. 

The stability threshold Q is the value of Q above which 
the disc is stable at all scales: 



STABLE Vfc 



g >■ 



(16) 



Q can be determined by imposing that the minimum of 
w^(fe) vanishes. Even though the calculations are very 
lengthy, the formula is simple, especially if expressed in 
terms of the 'right' parameters: 



Q = 2 



1/(2A) 



(17) 



where £t = I/^t is the Toomre scale. Eq. (16), with Q 
specified by Eq. (17), is our stability criterion. It reduces 
to Toomre's criterion Q > 1 in the limiting case of a non- 
turbulent disc: A = l, B = 2 (a = h = Q). 

The most unstable scale, i min — l/Z^min, is given by the 
formula: 



'C(g/Q) 



2/(B-A) 



(18) 



where Q/ Q measures the stability level of the disc (like Q in 
Toomre's case), and £ is the value of £min/^T at the stability 
threshold (Q/Q = 1): 



B-A 



A-i-\ iM 



(19) 



Eq. (18) reduces to £min/^T = | in Toomre's case. 

The growth rate of the most unstable scale, 7min = 
(— tjj^in)^^'^, is given by the formula: 



1 



(Q/Q) 



-2A/(B-A) 



(20) 



which vanishes at the stability threshold. Eq. (20) reduces 
to w^in/K^ = 1 — in Toomre's case. 



2.8 Case o > 1, 6 > i (1 -|- o) 

Even B < A < is a, Toomre-like case. In fact, although the 
scales at which pressure and rotation dominate are reversed, 
self-gravity still controls intermediate scales and oj'^{k) has 
a minimum, which determines whether the disc is stable or 
not. This means that, even now, the stability criterion is: 



STABLE Vfc 



Q>Q. 



(21) 



The stability threshold Q, the most unstable scale and its 
growth rate, fmin and 7inin, are given by the same formulae 
as in Eqs (17)-(20). 



The results of Sect. 2 are summarized in Fig. 1. The {a,b) 
plane is divided into four regions, where stability a la 
Toomre alternates with instability at small or large scales. 
The two shaded lines that separate those regions, and the 
point at which those lines intersect, represent transitions 
between different stability phases. Thus the corresponding 
stability criteria are hybrid (see Sects 2.2-2.4). The points 
(a, 6) = (0,0) and {a,b) = (— 1, i) represent the limiting 
case of a non-turbulent disc and the case investigated by 
Elmegreen (1996). To the best of our knowledge, this was 
the only theoretical work devoted to the gravitational insta- 
bility of clumpy /turbulent discs before ours. Fig. 1 also il- 
lustrates the most interesting stability regimes populated by 
models, sirrmlations and observations of astrophysical turbu- 
lence. Such points are discussed throughout the rest of this 
section. Specific predictions about strongly and weakly su- 
personic turbulence are then made in Sect. 3.1. 

The mass-size scaling relation, M oc Ei?^ oc tells 
us the natural bounds of a. In fact, D = a + 2 is the fractal 
dimension of the mass distribution, which ranges from to 
3, so we have — 2 < a < 1. Note that the upper bound 
corresponds to the case in which the destabilizing effect of 
self-gravity is scale- independent, i.e. to the vertical shaded 
line introduced above. 

Does even the other shaded line have a twofold mean- 
ing? Yes, and an important one! If the stabilizing effect of 
pressure has the same scale-dependence £is the effect of self- 
gravity, b = i (1 -I- a), then <x £Y1 oc M/£, which is the 
virial scaling relation. GMCs are then expected to clump 
along that line, i.e. to populate the transition regime be- 
tween stability a la Toomre and instability at small scales. 
For example, the well-known scaling laws E = constant and 
a oc ^^/^ (Larson 1981; Solomon et al. 1987) correspond 
to the point (a, &) = (0.^). Both Galactic and extragalac- 
tic GMCs show non-negligible dispersion around that point, 
especially along the virial line, as can be inferred from Bo- 
latto et al. (2008) and Heyer et al. (2009). The case of per- 
turbed galactic environments seems different. Rosolowsky & 
Blitz (2005) investigated the physical properties of GMCs 
in M64 (NGC 4826), an interacting molecule-rich galaxy. 



and found E oc m" '^*"'^ and a 



, which means 



(a, 6) = (5^3^,1.0 ± 0.3). If such scaling relations apply 
to individual GMCs, as was originally suggested, then each 
cloud is far from being in simple virial equilibrium. Besides, 
since (a, b) is below the virial line and on the right of the 
a = 1 line, the II2 disc is unstable at large scales (in the 
sense specified in Sect. 2.6) and the fractal dimension is for- 
mally higher than three. Alternatively, one may argue that 
those scaling relations arise from the superposition of more 
GMCs, each one characterized by the standard scaling laws, 
but with proportionality factors varying significantly over 
the disc (Rosolowsky, private communication). 

Now what about Hi? A turbulence model that is be- 
coming more and more popular is the one introduced by 
Fleck (1996), which predicts p^^^a oc f^^. To understand 
the meaning of this scaling relation, compare it with Kol- 
mogorov's law a oc f^^'^. Fleck's relation tells us that in a 
turbulent medium with both velocity and density fluctua- 
tions the quantity p^^^a plays a role similar to a in incom- 
pressible turbulence. Fleck (1996) assumed that S ~ pi. 
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Figure 1. The stability map of turbulence. The coordinates a and b are the exponents of the density- and velocity-size scaling relations 
[see Eqs (3) and (4)] . The shaded lines, their intersection point and the regions between them represent the variety of stability regimes 
possessed by clumpy /turbulent discs. The points (a, 6) = (0,0) and (a, 6) = (—1, ^) correspond to the limiting case of a non-turbulent 
disc and to the case investigated by Elmegreen (1996). The non-shaded part of the plane shows the natural range of a. The thick lines are 
phenomenological models of GMC and Hi turbulence, the semiclosed-head arrows represent observations, while the closed-head arrow 
represents simulations with increasing jVfach number. The solid circles correspond to points of special interest in astrophysical turbulence. 
Also shown are the cases illustrated in Figs 2 and 3 (hollow circles), and the degeneracy condition discussed in Sect. 4 (thin line). Note 
that the stability criterion is Q > 1 along such a line. 



which means I'^h (see Sect. 2.1). So his prediction cor- 
responds to the line b = | (2 — a), which crosses several 
stability regimes. The limiting case of Kolmogorov turbu- 
lence, {a,b) = (1, |), lies in the transition regime between 
stability a la Toomre and instability at large scales. Both 
high-resolution simulations of supersonic turbulence and H I 
observations populate the Toomre-like domain. Such simu- 
lations cluster along the Fleck line, near (a, 6) = (i, i), the 
case of Burgers turbulence (Kowal & Lazarian 2007; Krit- 
suk et al. 2007; Schmidt et al. 2008; Price & Federrath 2010; 
see also Fleck 1996 and references therein). In weakly super- 
sonic regimes, simulations cluster closer to the Kolmogorov 
limit (a, 6) = (1, i), as we will show in Sect. 3.1. Observed 
Hi intensity fluctuations, which are primarily due to co/c|3 
Hi (Lazarian & Pogosyan 2000), show large scatter around 
(a, 6) = (i, i), i.e. suggest a Kolmogorov scaling for both ve- 
locity and density fluctuations (e.g., Lazarian & Pogosyan 
2000; Elmegreen et al. 2001; Begum et al. 2006; Dutta et al. 
2008, 2009a, b). Such a scaling is also consistent with other 
Hi observations (e.g., Kim et al. 2007; Roy et al. 2008). The 



^ Hereafter we will omit 'cold' when referring to H I, since our 
analysis also focuses on the cold ISM (see Sect. 2.1). 



simulations by Agertz et al. (2009a), designed to explore 
the development of Hi turbulence in disc galaxies, suggest 
power-law indices consistent with the observed ones, except 
before the fragmentation of the disc (Agertz, Romeo et al., 
in preparation). 

Note that there is a very interesting and unexpected 
case where Fleck's model fits the observations well: the tran- 
sition from Hi to GMC turbulence. Clumps that climb up 
the Fleck line become progressively more self-gravitating 
and will virialize at {a,b) = (i, |). This is close to (a, 6) — 
(0, I), the reference point for GMC observations. The pre- 
dicted energy spectra are Es{k) oc dE'^/dfc « k~'^^^ and 
Ea{k) oc dcr^/dA: oc k~^^^^ (the Kolmogorov spectrum scales 
as fc-5/3). 

Low-resolution simulations of supersonic turbulence 
suggest an alternative scaling relation: p^^'^a oc (^^'^ (Kit- 
sionas et al. 2009). In comparison with the Fleck case, this 
scaling gives more weight to density fiuctuations and trans- 
lates into a steeper line: b = 5 (| — a), where it is again 
assumed that E ~ pi. The cases of Kolmogorov and Burgers 
turbulence correspond to (a, 6) = (1, |) and (a, 6) = (|, i), 
while the virialization point is (a, 6) = (^, |). The Kitsionas 
line crosses the same stability regimes as the Fleck line, but 
lies farther away from Hi and GMC observations. 
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3.1 Strongly vs. weakly supersonic turbulence 

Although our analysis focuses on strongly supersonic turbu- 
lence (see Sect. 2.1), here we extend it to weakly supersonic 
regimes (the case of a transonic or subsonic medium was 
considered in Sect. 2.1). 

How does the Mach number affect the density- and 
velocity-size scaling relations? And how does it affect the 
stability of the disc? To answer these questions, one should 
not compute a and b directly from the density and velocity 
spectral indices. One should first evaluate the typical density 
contrast of the medium (see Sect. 2.1). The density proba- 
bility distribution is approximately lognormal, and its mean 
/i and standard deviation SD depend on the rms Mach num- 
ber M (Padoan et al. 1997): fi = iSD^ ^ i ln(l + jM^). 
For such a distribution, the mass- weighted median density is 
Pmcd = pe^ ~ p{l + jM^Y^^ , where p is the mean density 
(see sect. 2.1.4 of McKee & Ostriker 2007). This provides 
a robust estimate of the typical density (mean plus fluctu- 
ations) in the medium. The corresponding density contrast 

is (5mod = (Pmcd - P)/P- 

In weakly supersonic turbulence [M ~ 2), we have 
5med ~ 0.4 so the mass column density at scale I is dom- 
inated by the mean density: Scfi ~ pi (hence a « 1) for 
^</i, while Eeff « ph (hence a « 0) for i^h. In contrast, 
the ID velocity dispersion at scale I is dominated by the tur- 

h, we can relate 6 to 
, and get (a, 6) ^ (1, i). 
Thus weakly supersonic 3D turbulence drives the disc to the 
boundary of the Toomre-like domain, near the Kolmogorov 
point. 

In strongly supersonic turbulence {A4 ^ 5) , we have 
^med ^ 2 so both the mass column density and the ID ve- 
locity dispersion are dominated by turbulent fluctuations. 
We can then compute a and b from the density and velocity 
spectral indices (see Sect. 2.1). Simulations show that in an 
isothermal medium with negligible self-gravity, with or with- 
out magnetic fields, the density spectrum flattens (Beres- 
nyak et al. 2005; Kim & Ryu 2005; Kowal et al. 2007) and the 
velocity spectrum steepens (e.g., Kowal & Lazarian 2007; 
Kritsuk et al. 2007; Price & Federrath 2010) as the Mach 
number increases. At Mach 7, we find (a, b) ~ (0.8, 0.4) from 
Kowal et al. (2007) and (a, b) ~ (0.3, 0.6) from Kowal & 
Lazarian (2007), using Fleck's model in both caseslfjln spite 
of the large uncertainties, it is clear that these points are 
on the left of (1, i) and, on average, closer to (i, i). Thus 
strongly supersonic 3D turbulence drives the disc well inside 
the Toomre-like domain, near the Burgers point. 



bulent term: a ?3 atur (£) (X r. For 
a using Fleck's model, b = ^ (2 — a] 



^ In the first case, dlni?p/dlnfc (i.e. — r) is given in table 2 of 
Kowal et al. (2007); a is computed from r as in Sect. 2.1 (^ </i): 
a = i (r + 1); and b is computed from a using Fleck's model: 
fe = i (2 - a). This yields (a, fe) = (0.75 ± 0.05,0.42 ± 0.02) for 
the sub-Alfvenic simulation and (a, b) = (0.80 ± 0.10, 0.40 ± 0.03) 
for the super-Alfvenic simulation. In the second case (of the main 
text), d In E„/d In k (i.e. — s) is given in fig. 1 of Kowal &; Lazarian 
(2007); b is computed from s as in Sect. 2.1: b = ^ (s — 1); and a 
is computed from b using Fleck's model: a = 2 — 3b. This yields 
(a, b) = (0.29 ±0.08, 0.57 ±0.03). This simulation is sub-Alfvenic, 
like the first one above. 
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Figure 2. The stability threshold of clumpy /turbulent discs, 
Q = Q{a,b,£o), where a and b are the exponents of the density- 
and velocity-size scaling relations, and £q is the typical scale at 
which density and velocity saturate. In addition, £t is the Toomre 
scale. The dashed lines and the hollow circles show the power-law 
behaviour predicted by Eq. (17) when density and velocity do not 
saturate; £o is then the fiducial scale at which those quantities are 
observed. The limiting case of a non-turbulent disc is Q = 1. 



4 ILLUSTRATIVE CASES 

Let us now illustrate the stability characteristics of 
clumpy/turbulent discs in a number of cases, those marked 
with hollow circles in the Toomre-like domain of our map 
(see Fig. 1). The points (i, i), (|, i) and (|, i) are typical 
values of (a, b) inferred from H I observations, high- and low- 
resolution simulations of supersonic turbulence. The points 
(|, 0) and (0, |) are the contributions of density and velocity 
fluctuations to the observed (a, b) = (^,^). For each case, 
we proceed in two ways. First, we assume that the density- 
and velocity-size scaling relations are perfect power laws, as 
given by Eqs (3) and (4), so that £o — 1/ko is the fiducial 
scale at which density and velocity are observed. We then 
evaluate the stability characteristics analytically using Eqs 
(17)-(20). Second, we consider more realistic scaling rela- 
tions, which take into account the saturation of density and 
velocity at large scales: 



1 



if ^ < £o , 
else ; 



1 else ; 



(22) 



(23) 



where £o is now the typical saturation scale. We then eval- 
uate the stability characteristics numerically using the dis- 
persion relation, Eq. (2), which we rewrite as: 



1 



+ 



V 



o! 

4 (£/£t) 



(24) 



where Q = noo/nGTia is the Toomre parameter and 
the Toomre scale. 
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Fig. 2 shows the stabihty threshold Q, i.e. the value of 
Q above which the disc is stable at all scales. The first curi- 
ous result is that such a diagnostic is highly degenerate. For 
example, look at the cases (a, 6) = (i, i) and (a, 6) = (|, i), 
which represent Hi observations and high-resolution simu- 
lations of supersonic turbulence. They have Q = 1\ Why do 
such cases degenerate into Toomre's case? Why doesn't tur- 
bulence show up? Eq. (17) gives us the answer: because the 
effects of density and velocity fluctuations cancel out when 
2 A — B = 0, i.e. along the line b = a (see the map). The 
cases (a, 6) = (|,0) and (a, fo) = (0, |) allow us to disen- 
tangle such effects. Density fluctuations that saturate at a 
typical scale £o tend to stabilize the disc by decreasing the 
stability threshold: Q < 1 if i'o > |^t, and Q = 1 otherwise. 
To understand this result, remember that such fluctuations 
reduce the self-gravity term in the dispersion relation by a 
factor [see Eqs (22) and (24)], and that self-gravity is 
destabilizing. Velocity fluctuations have an antagonistic ef- 
fect. They reduce pressure by a factor V| [see Eqs (23) and 
(24)], and tend to destabilize the disc by increasing the sta- 
bility threshold: Q > 1 if ^o^|^t, and Q = 1 otherwise. 
When density /velocity fluctuations do not saturate, their 
effect becomes destabilizing/stabilizing if the fiducial scale 
is small. 

Fig. 3 shows the most unstable scale, £min, for two values 
of Q/Q. Remember that this quantity measures the stabil- 
ity level of the disc, like Q in Toomre's case. So Q/ Q = 1 
(top panel) means that the disc is marginally unstable, while 
Q/ Q = 0.7 (bottom panel) means that the disc is mod- 
erately unstable. In contrast to Q, Imin is not degenerate. 
Turbulence now has a significant effect in the case of H I 
observations, since the contributions of density and velocity 
fluctuations are no longer antagonistic. For £o ~ It, £ min is 
about 30 to 50 per cent smaller than in Toomre's case, de- 
pending on the value of Q/ Q. The impact of turbulence is 
stronger in the case of high-resolution simulations: the most 
unstable scale is typically half an order of magnitude below 
the expected value. Turbulent effects become less important 
at low resolution. 

The growth rate of the most unstable scale, 7min, is 
independent of £o and vanishes for Q/ Q = 1 [see Eq. (20)]. 
For Q/ Q < 1, the effects of density and velocity fluctuations 
cancel out when A/{B — A) = 1, i.e. b = a. This degeneracy 
condition is the same as that found for Q, and has the same 
consequences. 



J2o/-2t 

Figure 3. The most unstable scale of clumpy/turbulent discs, 
^min = ^min(ii f*! ^Oi Q/ Q) , where a and b are the exponents of 
the density- and velocity-size scaling relations, £o is the typical 
scale at which density and velocity saturate, and Q/Q is the 
stability level. In addition, £x is the Toomre scale. The dashed 
lines show the power-law behaviour predicted by Eqs (18) and 
(19) when density and velocity do not saturate; £q is then the 
fiducial scale at which those quantities are observed. Also shown 
is the limiting case of a non-turbulcnt disc (dotted line) . 



5 CONCLUSIONS 

• Observations and simulations of the interstellar medium 
are revealing its turbulent nature with higher and higher fi- 
delity. Such information must then be taken into account 
when analysing the stability of galactic discs. Our contri- 
bution is a natural extension to Toomre's work, which will 
prove useful for both low- and high-redshift analyses. 

• Turbulence has a deep impact on the gravitational in- 
stability of the disc. It excites a rich variety of stability 
regimes, several of which have no classical counterpart. We 
illustrate this result in the form of a map, which relates our 
stability scenario to the phenomenology of interstellar tur- 
bulence: GMC/Hl observations, simulations and models. 

• GMC turbulence drives the disc to a regime of tran- 
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sition between instability at small scales and stability a 
la Toomre. Toomre's criterion works instead typically well 
when applied to discs of cold Hi, since the effects of density 
and velocity fluctuations tend to cancel out. Even so, H I 
turbulence produces a clear signature in disc morphology. It 
reduces the characteristic scale of instability by 30-50 per 
cent or more, depending on the value of Q and on the shape 
of the energy spectrum. The transition from Hi to GMC 
turbulence occurs when E ~ £^^^ and a ~ ^^/^ (for more 
information see Sect. 3). 

• Coming astronomical facilities such as will be 

able to resolve the scaling properties of galactic turbulence 
up to very high redshifts. Using our map, such data will 
show up as evolutionary tracks, which will reveal the inter- 
play between gravitational instability and turbulence during 
the galaxy life. In turn, this will be useful for constraining 
the sources of galactic turbulence and for understanding the 
evolution of cosmic star formation. 
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